Bursts in a fiber bundle model with continuous damage 
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We study the constitutive behaviour, the damage process, and the properties of bursts in the 
continuous damage fiber bundle model introduced recently. Depending on its two parameters, 
the model provides various types of constitutive behaviours including also macroscopic plasticity. 
Analytic results are obtained to characterize the damage process along the plastic plateau under 
strain controlled loading, furthermore, for stress controlled experiments we develop a simulation 
technique and explore numerically the distribution of bursts of fiber breaks assuming infinite range 
of interaction. Simulations revealed that under certain conditions power law distribution of bursts 
arises with an exponent significantly different from the mean field exponent 5/2. A phase diagram 
of the model characterizing the possible burst distributions is constructed. 

PACS number(s): 46.50.+a, 62.20.Fe, 62.20.Mk 
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I. INTRODUCTION 



Recently, the breakdown of disordered materials under 
externally imposed stresses has attracted much attention 
and by now several aspects of the breakdown process are 
well understood ||l|,g| . The possibilities of pure analytical 
approaches for breakdown phenomena are rather limited, 
hence, computer simulation is an indispensable tool in 
this field. Models, computer simulations are based on, 
can be classified as lattice models and fiber bundle mod- 
els. In lattice models the elastic medium is represented 
by a spring (beam) network, and disorder is captured ei- 
ther by random dilution or by assigning random failure 
thresholds to the bonds M. The failure rule usually ap- 
plied in lattice models is discontinuous and irreversible: 
when the local load exceeds the failure threshold of a 
bond, the bond is removed from the calculations {i.e. its 
elastic modulus is set to zero). 

A very important class of models of material failure 
are the fiber bundle models (FBM) [^-Q, which have 
been extensively studied during the past years. These 
models consists of a set of parallel fibers having statisti- 
cally distributed strength. The sample is loaded parallel 
to the fibers direction, and the fibers fail if the load on 
them exceeds their threshold value. In stress controlled 
experiments, after each fiber failure the load carried by 
the broken fiber is redistributed among the intact ones. 
The behaviour of a fiber bundle under external loading 
strongly depends on the range of interaction, i. e. on the 
range of load sharing among fibers. Exact analytic re- 
sults on FBM have been achieved in the framework of 
the mean field approach, or global load sharing, which 
means that after each fiber breaking the stress is equally 



distributed on the intact fibers implying an infinite range 
of interaction and a neglect of stress enhancement in the 
vicinity of failed regions p-0,Eq-p8[. In spite of their 
simplicity FBM's capture the most important aspects of 
material damage and they provide a deep insight into the 
fracture process. Over the past years several extensions 
of FBM have been carried out by considering stress local- 
ization (local load transfer) ||9|, p^ , p^,p^ , 32 p4|, the effect 
of matrix material between fibers |1^ 17,|29 EM], possi- 
ble non- linear behavior of fibers g , thermally activated 
breakdown [Q of fibers, and coupling to an elastic block 

il- 

Very recently, a novel continuous damage law has been 
introduced in lattice models Pq] of fracture. In this 
model when the failure threshold of a lattice bond is ex- 
ceeded the elastic modulus of the bond is reduced by 
a factor a (0 < a < 1), furthermore, multiple failures 
of bonds are allowed. Simulations revealed that under 
strain controlled loading the system develops into a self 
organized state, which is macroscopically plastic, and is 
characterized by a power law distribution of avalanches 
of breaks. We worked out an extensions of fiber bundle 
models by implementing a continuous damage law for the 
fibers ||l^ , in the spirit of Ref. |^ . It has been demon- 
strated in Ref. 113] that the continuous damage fiber bun- 
dle model (CDFBM) provides a broad spectrum of de- 
scription of materials varying its parameters and for cer- 
tain parameter settings the model recovers the variants 
of FBM known in the literature. CDFBM can be relevant 
for materials where the microscopic damage mechanism 
is gradual multiple failure of components, i.e. matrix 
and fibers |39 41]. Very recently, the CDFBM was fur- 
ther developed by Moral et al. taking into account time 
dependence in the failure process |l3]. 
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One of the most appealing results on CDFBM was that 
the multiple failure of brittle elements can give rise to a 
macroscopic plastic behaviour of the specimen, which is 
then followed by a hardening or softening regime, fur- 
thermore, under certain conditions damage localization 
occurs. However, the microscopic damage process of 
CDFBM has not been explored. The main goal of the 
present paper is to reveal the microscopic failure pro- 
cess in order to understand the emergence of the plas- 
tic macroscopic state. Analytic results are obtained to 
characterize the damage process along the plateau un- 
der strain controlled loading, furthermore, for stress con- 
trolled experiments we develop a simulation technique 
and explore numerically the distribution of bursts of fiber 
breaks. The effect of localization on the process of dam- 
age is clarified. A phase diagram of the model charac- 
terizing the possible constitutive behaviours and burst 
distributions is constructed in terms of the two parame- 
ters of the model. 



II. MODEL 



The continuous damage fiber bundle model is an ex- 
tension of the commonly used fiber bundle models by 
generalizing the damage law of fibers. The model system 
is composed of A^ parallel fibers with identical Young- 
modulus Ef but with random failure thresholds di, i = 
1, . . . , A. The failure strength di of individual fibers is 
an independent identically distributed random variable 
with a probability density p{d) and a cumulative prob- 
ability distribution P{d) ~ L p{x)dx. The fibers are 
assumed to have linear elastic behaviour up to breaking 
(brittle failure) . Under uniaxial loading of the specimen 
a fiber fails if it experiences a load larger than its break- 
ing threshold di. In the framework of our model at the 
failure point the stiffness of the fiber gets reduced by a 
factor a, where < a < 1, i.e. the stiffness of the fiber 
after failure is aEf. In principle, a fiber can now fail 
more than once and the maximum number kmax of fail- 
ures allowed for fibers is a parameter of the model. Once 
a fiber has failed its damage threshold di can either be 
kept constant for the further breakings (quenched dis- 
order) or new failure thresholds of the same distribution 
can be chosen (annealed disorder), which can model some 
microscopic rearrangement of the material after failure. 
The damage law of the model is illustrated in Fig. |] for 
both types of disorder. The characterization of dam- 
age by a continuous parameter corresponds to describe 
the system on length scales larger than the typical crack 
size. This can be interpreted such that the smallest ele- 
ments of the model are fibers and the continuous damage 
is due to cracking inside fibers. However, the model can 
also be considered as the discretization of the system on 
length scales larger than the size of single fibers, so that 
one element of the model consists of a collection of fibers 
with matrix material in between. In this case the micro- 



scopic damage mechanism resulting in multiple failure of 
the elements is the gradual cracking of matrix and the 
breaking of fibers. In the following we refer to the ele- 
ments of the continuous damage FBM as fibers, but we 
have the above two possible interpretations in mind. 
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FIG. 1. The damage law of a single fiber of the contin- 
uous damage model when multiple failure is allowed a) for 
quenched, and b) for annealed disorder. The horizontal lines 
indicate the damage threshold di. 



After failure the fiber skips a certain amount of load 
which has to be taken by the other fibers. For the 
load redistribution we assume infinite range of interaction 
among fibers (mean field approach), furthermore, equal 
strain condition is imposed which implies that stiffer 
fibers of the system carry more load. At a strain e the 
load of fiber i that has failed k{i) times reads as 



Me)=Efa'^'>e 



(1) 



where Eja^^^^ is the actual stiffness of fiber i. It is impor- 
tant to note that, in spite of the infinite interaction range, 
Eq. (||) is different from the usual global load sharing 
where all the intact fibers carry always the same amount 
of load. In the following the initial fiber stiffness Ef will 
be set to unity. 



III. CONSTITUTIVE LAWS 



Here we provide a derivation of the constitutive law 
of the continuous damage FBM in a more transparent 



way than in Ref. [|2J. This general theoretical frame- 
work facilitates to obtain analytic results also for the mi- 
croscopic failure process. The key quantity is the prob- 
ability pbk (e) that during the loading of a specimen an 
arbitrarily chosen fiber failed precisely /c-times at a strain 
e, where k = 0, . . . , kmax denotes the failure index, and 
A: = is assigned to the intact fibers. pbk{e) can be cast 
in the following form for annealed disorder 



fe-i 



pbk{e)=[l-P{a''e)]l[P{a^e), 



(2) 



for < A: < fc; 



max 

krr 



1, (3) 



and p&fe„_(e)= [] ^(«'^)' 

and for quenched disorder 

pbois) = 1 - Pis), 

pbkie) = P{a^^'^e) - Pia'^e), for 1 < k < k„ 

and pfefc„„Je) = P(a'''""--ie). 



It can be easily seen that the probabilities Eqs. (|2|J^) 
fulfill the normalization condition ^^."q'" p6fe(e) — 1. Av- 
erage quantities of the fiber ensemble during a loading 
process can be calculated using the probabilities Eqs. 
(^,||). For instance, the average load on a fiber F/N 
at a given strain e reads as 



F 

N 



J2 '^^P^kie) 



fc=0 



(4) 



which provides the macroscopic constitutive behaviour 
of the model, and the expression in the brackets can be 
considered as the macroscopic effective Young modulus 
of the sample {Ef = 1). The single terms in the sum give 
the load carried by the subset of fibers of failure index k. 
The variants of fiber bundle models used widespread in 
the literature can be recovered by special choices of the 
parameters kmax and a of the model. A micromechanical 
model of composites ||l(^ , |l7| , po| , pi| can be obtained with 
the parameter values kmax = 1, a 7^ 



F 

N 



£ [1 - Pie)] + aeP{e), 



(5) 



while setting kmax — 1, a — 0, i.e. skipping the second 
term in Eq. (1^) results in the classical dry bundle model 
of Daniels j||] . 

In Fig. ^we show the explicit form of the constitutive 
law with annealed disorder for different values of kmax in 
the case of the WeibuU distribution 



P(d) = l-exp(-(d/4)'"), 



(6) 



where m is the WeibuU modulus and dc denotes the char- 
acteristic strength of fibers. The parameter values are set 
to m = 2, dc ~ 1 in all the calculations. 



Note that in the constitutive equation Eq. (Q) the term 
of the highest failure index k^ax can be conceived such 
that the fibers have a residual stiffness of a'"'™"^ after hav- 
ing failed kmax times. This residual stiffness results in a 
hardening of the material, hence, the F/N curves in Fig. 
11(a) asymptotically tend to straight lines with a slope 
gkman Increasing kmax the hardening part of the consti- 
tutive behavior is preceded by a longer and longer plastic 
plateau, and in the limiting case of kmax —^ c)o the ma- 
terials behavior becomes completely plastic. A similar 
plateau and asymptotic linear hardening has been ob- 
served in brittle matrix composites, where the multiple 
cracking of matrix turned to be responsible for the rela- 
tively broad plateau of the constitutive behavior, and the 
asymptotic linear part is due to the linear elastic behav- 
ior of fibers remained intact after matrix cracking [|39| . 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

e 

FIG. 2. Constitutive behaviour of the model of annealed 
disorder a) with b) without residual stiffness at a = 0.8 for 
different values of kmax- In b) the lowest curve presents the 
constitutive behaviour of the dry bundle model for compari- 
son. 



In order to describe macroscopic cracking and global 
failure of a specimen instead of hardening, the resid- 
ual stiffness of the fibers has to be set to zero after a 
maximum number k* 



of allowed failures |12 35[|. In this 
case the constitutive law can be obtained from the gen- 
eral form Eq. (m by replacing kmax in the upper limit 



of the sum by fc* — 1. A comparison of the constitu- 
tive laws of the dry and continuous damage FBM with 
global failure is presented in Fig. |2|(b). One can observe 
that the dry FBM constitutive law has a relatively sharp 
maximum, while the continuous damage FBM curves ex- 
hibit a plateau whose length increases with increasing 
k* . Note that the maximum value of F/N corresponds 
to the macroscopic strength of the material, furthermore, 
in stress controlled experiments the plateau and the de- 
creasing part of the curves cannot be reached. However, 
by controlling the strain e, the plateau and the decreas- 
ing regime can also be realized. The value of the driving 
stress a = F/N corresponding to the plastic plateau, and 
the length of the plateau are determined by the damage 
parameter a, and by kmax, k*: Decreasing a at a fixed 
kmax, k* , or increasing kmax, k* at a fixed a give rise to 
an increase of the plateau's length. 

In Fig. y we compare the constitutive laws of the model 
with different types of disorder for hardening and soften- 
ing of fibers. It can be seen that there is no qualitative 
difference between the curves of annealed and quenched 
disorder, however, the inset shows that when the fibers do 
not have residual stiffness (softening) the local shape of 
the curves around the maximum is different. The differ- 
ent types of maxima of the constitutive curve of the two 
different disorders result in a very interesting behaviour 
of the burst distributions, which will be discussed later. 
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FIG. 3. Comparison of the constitutive behaviours with 
annealed and quenched disorder. We show data with and 
without remaining stiffness. The inset demonstrates how the 
shape of the constitutive curves changes when increasing k* 
with different types of disorders. 



IV. DAMAGE 



The damage state of the model at a certain e can be 
characterized by the average number of failures occurred. 



Based on the probabilities Eqs. 
damage variable D{e) as 



(g,^), we introduce a 



Die) 



1 



k ^-^ 

fc— 1 



kpbk{e) 



(7) 



which is an integral quantity of the damage process. 
From the properties of pbk{e) it can be seen that D 
is a monotonically increasing function, and D G [0,1]. 
Then the average number of failures can be obtained as 
NkmaxD{e). Fig. ^ illustrates the behaviour of D for 
three different values oi kmnx- 




FIG. 4. The damage variable D for annealed (open sym- 
bols) and quenched (filled symbols) disorder for several dif- 
ferent values of kmax- The damage variable was chosen to be 
a = 0.8. 



It can be observed in Fig. Q that the overall behaviour 
of the damage variable D is nearly the same for annealed 
and quenched disorder, however, there is a significant dif- 
ference between the microscopic damage processes in the 
two cases. In spite of the infinite range of interaction 
among fibers, localization of damage occurs for the case 
of quenched disorder. It means that weaker fibers tend to 
break more often than the stronger ones. For quenched 
disorder, the strain £„ where the weakest fiber of failure 
threshold d„i reaches kmax, is Sm = rfr?j./a'^"""°- Hence, 
at this loading stage the failure index k of fibers as a 
function of the damage threshold d can be obtained as 



k{d) 



— In — 
In a dm. 



(8) 



Localization of damage means that fc is a decreasing func- 
tion of d, and it can be seen from Eq. (|g) that the localiza- 
tion gets more pronounced when the damage parameter 
a ^ 1. This localization effect is illustrated in Fig. ||, 
where the analytic result Eq. (p[) is compared to simula- 
tions for three different values of a. 
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FIG. 5. The failure index k of fibers as a function of their 
failure threshold d for several different values of the damage 
parameter a. The number of fibers was chosen to be A'^ = 500. 
The continuous lines represent the corresponding analytical 
results of Eq. (S). 



we study the distribution of bursts in our CDFBM under 
strain and stress controlled conditions. 



A. Strain controlled case 



Under strain controlled conditions of fiber bundles 
there is no load transfer from broken to intact fibers since 
the load carried by each fiber is determined by the exter- 
nally imposed strain and the local fiber stiffness accord- 
ing to Eq. (nf). This implies that the number of fibers 
which break due to an infinitesimal increase of the ex- 
ternal strain is completely determined by the statistics 
of fiber strength, i.e. by p{d) and P{d). It has been 



discussed in Sec. ^ that the plastic plateau and the de- 
creasing part of the constitutive law can only be realized 
in strain controlled experiments. To reveal the nature 
of ductility arising in our model it turns to be useful to 
study the statistics of bursts occurring under strain con- 
trolled conditions. 



V. DISTRIBUTION OF BURSTS 

One of the most interesting aspects of the damage 
mechanism of disordered solids is that the breakdown is 
preceded by an intensive precursor activity in the form of 
avalanches of microscopic breaking events |g^-^,^,^ . 
Under a given external load F a certain fraction of fibers 
fails immediately. Due to the load transfer from bro- 
ken to intact fibers this primary fiber breaking may ini- 
tiate secondary breaking that may also trigger a whole 
avalanche of breakings. If F is large enough the avalanche 
does not stop and the material fails catastrophically. For 
the dry FBM it has been shown by analytic means that 
in the case of global load transfer the size distribution of 
avalanches follows asymptotically a universal power law 
with an exponent —5/2 p9| , |30| , p8{ , however, in the case of 
local load transfer no universal behaviour exists, and the 
avalanche characteristic size is bounded p^. This pre- 
cursory activity can also be observed experimentally by 
means of the acoustic emission analysis. Acoustic emis- 
sion measurements have revealed that for a broad variety 
of disordered materials the response to an increasing ex- 
ternal load takes place in bursts having power law size 
distribution over a wide range ||36|-|38|] . 

Introducing a continuous damage law in lattice models, 
simulations revealed that under strain controlled condi- 
tions the system tends to a steady state, which is macro- 
scopically plastic Q, similarly to our case. Due to the 
long range interaction, the plastic steady state is charac- 
terized by power law distributed avalanches of breaks and 
it has been argued that the underlying damage mecha- 
nism displays self organized criticality. In the following 
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FIG. 6. kmaxpbtot as a function of e for annealed disorder, 
comparison of simulations and analytic results of Eq. (Mj (con- 
tinuous lines) . The integral of the functions is always equal to 
k* . In the upper part of figure the corresponding constitutive 
curves are also presented for comparison. 



The basic quantity to characterize bursts is the proba- 
bihty pbi.'^^{e)de that a fiber, which has failed k times up 
to strain e imposed externally, will fail again under an in- 
finitesimal strain increment de. From Eqs. (||,y) p6^^^(e) 
can be cast in the form for annealed disorder 






^ ^ ^p{a'e)a' 



^^ Pia^e) ' 



(9) 



and for quenched disorder 

pbl+\e)^p{a''e)a'', k = Q,.. 



h — 1 



(10) 



and the total probability of fiber breaking can be ob- 
tained by summing over k 



pbtotie) 



1 



E Pbt^'i^)- 



(11) 



fe=o 




The number of fiber failures occurring in the strain 
interval [e, e -I- de] can be obtained as NkmaxPbtot{()de. 
This is a very important characteristic quantity of the 
microscopic damage process since it can be monitored 
experimentally by means of acoustic emission techniques. 
The behaviour oipbtot is shown in Fig. |6|for the softening 
case with several values of k* , where also the correspond- 
ing constitutive curves are presented. It can be seen that 
pbtot lias a maximum where the plastic regime of the con- 
stitutive curve starts, and it is a decreasing function of 
e in the whole plastic region. Due to the stiffness reduc- 
tion of the system caused by the subsequent failures, in 
the plastic regime the same increase of strain results in 
smaller and smaller load increments on fibers, and hence, 
pbtot and the number of failures decreases. It also implies 
that the breaking activity, which can be measured by 
acoustic emission techniques, decreases along the plateau 
in agreement with experiments pj-Esl . 

It follows from the above argument that decreasing 
the value of the damage parameter a while kmax is kept 
fixed, the length of the plastic plateau, preceding the de- 
creasing or hardening part of the constitutive behaviour, 
increases since larger strain is required to achieve succes- 
sive failure. This is demonstrated in Fig. ui where one 
can also see that for small a the constitutive curve devel- 
ops distinct maxima. In order to clarify the occurrence 
of these maxima in the plastic plateau, in Fig. M we also 
plotted pbil^^ for three different values of fc at a = 0.4. 
With decreasing a the length of the plastic plateau in- 
creases, however, the consecutive maxima of pb'^'^^ get 
more and more separated giving rise to visible maxima 
in the plateau. The broader the disorder distribution is, 
the smaller the value of a is where the maxima of F/N 
appear. 
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FIG. 7. a) The constitutive behaviour varying the damage 
threshold at a fix ifc*. fe) pfe^+^(e) for a = 0.4. 



FIG. 8. The energy dissipation rate for two diflerent values 



of k„ 



The energy dissipation rate is also a very important 
aspect of the ductile regime of the model. The energy 
dissipation rate Edis (e) is defined so that the energy dis- 
sipated due to the failure of fibers in the strain interval 
[e, e + de] can be obtained as Edis{^)d€ 



Edis{^) — 



E 

fc=0 



1 



e^a^l-a) 



Pb^'i^l 



(12) 



where the expression in the bracket provides the energy 
dissipated by the failure of a fiber which has already failed 
fc-times. In Fig. g the energy dissipation rate Edis{f) is 
plotted for two different values of kmax- Comparing Fig. 
H to the corresponding constitutive curves in Fig. it can 
be observed that in the plastic regime Edis{t) is constant. 



B. Stress controlled case 

Under stress controlled loading conditions the micro- 
scopic dynamics of the damage process is more compli- 
cated than in the strain controlled case, since the fail- 
ure of each fiber is followed by a redistribution of load, 
which can provoke further fiber breakings resulting in 
an avalanche of failure events. Studying the statistics of 
avalanches under quasi-static loading of a specimen, im- 
portant information can be gained about the dynamics 
of damage, which can be then compared to the results 
of acoustic emission experiments. Due to the difficulties 
of the analytic treatment, we develop a simulation tech- 
nique and explore numerically the properties of bursts in 
our continuous damage fiber bundle model. The inter- 
action of fibers, the way of load redistribution is crucial 
for the avalanche activity. A very important property 
of CDFBM is that in spite of the infinite range of inter- 
action the load on intact fibers is not equal, but stiffer 
fibers cary more load, furthermore, for quenched disorder 
damage localization occurs, which might affect also the 
avalanche activity. 

To implement the quasi-static loading of a specimen of 
TV fibers in the framework of CDFBM, the local load on 
the fibers fi has to be expressed in terms of the external 
driving F. Making use of Eq. (p it follows that 



N 



N 



p^Ef^-^j:-' 



(0 



(13) 



and hence, the strain and the local load on fibers can be 
obtained as 



F 



Eti «''^« ^ 



f^^F 
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Ef=i «'=« 



(14) 



when the external load F is controlled. The simulation 
of the quasi-static loading proceeds as follows: in a given 
stable state of the system we determine the load on the 
fibers fi from the external load F using Eq. (^#. The 
next fiber to break can be found as 



di 

TV 



r > 1, 



(15) 



i.e. that fiber breaks for which the ratio di/ fi is the 
smallest. Here i* denotes the index of the fiber to break, 
di is the damage threshold of fiber i, and fi is the lo- 
cal load on it. To ensure that the local load of a fiber 
is proportional to its stiffness, the external load has to 
be increased in a multiplicative way, so that F — > rF is 
imposed, and the failure index of fiber i* is increased by 
one k{i*) — > k{i*) + 1. After the breaking of fiber i*, 
the load fi carried by the fibers has to be recalculated 
making use of Eq. ([l4|), which provides also the correct 
load redistribution of the model. If there are fibers in 
the state obtained, whose load exceeds the local breaking 
threshold, they fail, i.e. their failure index is increased 
by 1 and the local load is again recalculated until a stable 
state is obtained. A fiber cannot break any longer if its 
failure index k has reached k* or k^ax during the course 
of the simulations. This dynamics gives rise to a complex 
avalanche activity of fiber breaks, which is also affected 
by the type of disorder. The size of an avalanche 5* is 
defined as the number of breakings initiated by a single 
failure due to an external load increment. 

Simulations revealed that varying the two parameters 
of the model kmax, a, or fc*, a and the type of disorder, the 
CDFBM shows an interesting variety of avalanche activ- 
ities, characterized by different shapes of the avalanche 
size distributions. 
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FIG. 9. Avalanche size distributions for different values of 
kmax and a when fibers have remaining stiffness and the disor- 
der is annealed. The number of fibers was A*' = 1600 and aver- 
ages were made over 2000 samples. The number of avalanches 
n of size 5* are shown to demonstrate also how the total num- 
ber of avalanches changes. 



In Fig. g the histograms n{S) of the avalanche sizes S 
are shown which were obtained for a system of remaining 
stiffness and annealed disorder with Weibull parameters 
m — 2, dc = 1. Since in the limiting case of a — > the 
CDFBM recovers the global load sharing dry fiber bun- 
dle model, in Fig. ^ the curves with small a and k^ax = 1 
are power laws with an exponent a = 5/2 in agreement 
with the analytic results |g^,Q. Increasing the value of 
a at a fixed kmax only gives rise to a larger number of 
avalanches, i.e. parallel straight lines are obtained on a 
double logarithmic plot, but the functional form of ^(5") 
does not change. However, when a exceeds a critical 
value Oc (cc ~ 0.3 was obtained with the Weibull param- 
eters specified above) the avalanche statistics drastically 
changes. At a fixed a > ac when kmax is smaller than a 
specific value kc{a), the avalanche sizes show exponential 
distribution, while above kc{a) the distribution takes a 
power law form with an exponent f3 = 2.12 ± 0.05. 



law distribution of avalanches with a constant exponent 
/3 = 2. 12 ±0.05 significantly different from the mean field 
exponent a — 5/2 |29|3^. It is important to emphasize 
that the overall shape of phase diagram is independent of 
the type of the disorder (annealed or quenched), more- 
over, the specific values Oc ~ 0.3 and kc{a) depend on 
the details of the disorder distribution p{d). 

A very different behavior was obtained for the system 
when fibers do not have remaining stiffnes after a k* num- 
ber of failures. Simulations revealed that in this case the 
avalanche statistics strongly depends on the type of disor- 
der. When the disorder is quenched the size distribution 
of avalanches follows always the dry bundle results for 
the whole domain of parameters, i.e. N(S) shows power 
law behaviour with an exponent a = 5/2. When A;* > 1 
the larger number of breakings results in more avalanches 
but the overall distribution does not change. 
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FIG. 10. Phase diagram for the continuous damage model 
with remaining stiffness for both types of disorder. The func- 
tional form of the avalanche statistics is given in the param- 
eter regimes. The location of the Dry Bundle Model (DBM) 
in the parameter space is also indicated. 



Based on the above results of simulations a phase 
diagram is constructed which summarizes the proper- 
ties of avalanches with respect to the parameters of the 
model. Fig. |l^ demonstrates the existence of three dif- 
ferent regimes. If the damage parameter a is smaller 
than Oc, the dynamics of avalanches is close to the sim- 
ple Dry Bundle Model characterized by a power law of 
the mean field exponent a = 5/2. However, for a > Qc 
the avalanche size distribution depends on the number 
of failures kmax alowed. The curve of kc{a) in the phase 
diagram separates two different regimes. For the parame- 
ter regime below the curve, avalanche distributions with 
an exponential shape were obtained. However, the pa- 
rameter regime above kc{a) is characterized by a power 
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FIG. 11. Avalanche size distributions for different values of 
k* at a fixed a > ac when fibers have no remaining stiffness 
and the disorder is annelaed. 



Nevertheless, when the disorder is annealed the system 
shows a more complex behaviour. When a falls below a 
certain critical value Oc the results are similar to DBM 
independently of the value of fc*, however, for a > Oc a 
novel avalanche dynamics appears (for the present values 
of the Weibull parameters Qc ~ 0.35 was obtained). In 
Fig. [ll| the avalanche distributions are shown for an a 
value above ac, varying the value of k*. It is very im- 
portant to emphasize that the curves in all the cases can 
be well fitted with a power law, however, the value of 
the exponent depends on fc*. Two extreme cases can be 
distingushed: for fc* = 1 the system recovers the DBM 
avalanche dynamics. On the other hand, for fc* > kc{a) 
the exponent of the power laws is (3 = 2.12 ± 0.05, simi- 
larly to the case of remaining stiffness. Below kc{a) the 



exponents vary as a function of k* between the mean field 
exponent a and /3. 

Based on Refs. [ p9pO| ] the different types of avalanche 
size distributions can also be understood up to some 
extent in terms of the constitutive curves of Sec. 111. 
Comparing Fig. ^ and Figs. ^ ^ it can be recognized 
that if the constitutive curve has a single quadratic max- 
imum the corresponding avalanche size distribution of 
CDFBM follows the mean field results, while other types 
of avalanche statistics arises when this condition does not 
hold. 
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VI. CONCLUSIONS 



A detailed analytical and numerical study of the con- 
tinuous damage fiber bundle model is presented. The 
model is an extension of the classical fiber bundle model 
by introducing a continuous damage law, and allowing 
for multiple failure of fibers with quenched and annealed 
disorders. A simple general derivation of the constitu- 
tive behaviour of the model is provided, which also fa- 
cilitates to obtain analytic results for the microscopic 
damage process. Varying its parameters, the model pro- 
vides a broad spectrum of description of materials rang- 
ing from strain hardening to perfect plasticity, and hence, 
the model can be relevant to describe the damage process 

43[ |. It is a remarkable 
e failure of brittle ele- 



tip 



of various types of materials | 

feature of the model that mu 

ments can result in a macroscopically plastic state, which 

has also been observed experimentally in materials where 

the damage mechanism is the gradual multiple failure of 

ingredients [|2|J4|]. 

The present study focused on the microscopic dam- 
age process to understand the emergence of the plastic 
plateau under strain controlled loading, and the resulted 
avalanche activity under stress controlled loading of the 
continuous damage fiber bundle model. Analytic results 
are obtained to characterize the damage process along 
the plateau under strain controlled loading, furthermore, 
for stress controlled experiments a simulation technique 
was developed and the distribution of avalanches of fiber 
breaks was explored numerically. Simulations showed 
that depending on the parameters of the model the dis- 
tribution of bursts of fiber breaks can be exponential or 
power law. Based on extensive computer simulations, a 
phase diagram characterizing the possible avalanche dis- 
tributions is constructed in terms of the two parameters 
of the model. One of the most appealing outcomes is 
that the model has a broad parameter regime where the 
avalanche statistics shows a power law behaviour with an 
exponent significantly different from the well know mean 
field exponent, in spite of the infinite range of interaction 
among fibers. The results obtained have relevance to un- 
derstand the acoustic emission measurements performed 
on various elasto-plastic materials 
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